Polaron formation in the presence of Rashba spin-orbit coupling: implications for 

spintronics 
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We study the effects of the Rashba spin-orbit coupling on the polaron formation, using a suitable 
generalization of the Momentum Average approximation. Contrary to previous investigations of this 
problem, we find that the spin-orbit interaction decreases the effective electron-phonon coupling. It 
is thus possible to lower the effective mass of the polaron by increasing the spin-orbit coupling. We 
also show that when the spin-orbit coupling is large as compared to the phonon energy, the polaron 
retains only one of the spin polarized bands in its coherent spectrum. This has major implications 
for the propagation of spin-polarized currents in such materials, and thus for spintronic applications. 



PACS numbers: 71.38.-k, 71.70.Ej 

In the field of spintronics there is a continued effort 
to develop means of efficiently manipulating the spin of 
electrons [l| . One widely studied approach is to use spin- 
orbit coupling (SO) to lift the spin degeneracy, for exam- 
ple the Rashba SO coupling Q that arises in a confined 
system if the quantum well lacks inversion symmetry. Ex- 
perimentally, the Rashba effect has been seen in many 
systems, 6.17. semiconductor heterostructures like GaAs 
and InAs, surface states of metals like Au(lll)[3| and 
surface alloys Hkc Bi/Ag [3| or Pb/Ag 

Confined two-dimensional (2D) systems may also cou- 
ple strongly to optical phonons of the substrate. Tuning 
of these electron-phonon (el-ph) interactions was shown 
to be experimentally viable in organic single-crystal tran- 
sistors [6| . Strong el-ph coupling is interesting because it 
leads to polaron creation, whereas the coherent quasipar- 
ticle is an electron surrounded by a phonon cloud. The 
polaron dispersion and effective mass can be significantly 
renormalized from those of the bare band electron Q]. 

An interesting question is whether the el-ph and the 
SO coupling can be used in conjuction to tailor differently 
the properties of the two bands with different spin. The 
interplay of the Rashba SO and the el-ph interactions has 
been considered previously for continuous systems with 
parabolic bands and weak el-ph coupling, using the self- 
consistent Born approximation Q . The main conclusion 
was that SO enhances the effective el-ph coupling due to 
a topological modification of the Fermi surface, namely 
an effective reduced dimensionality of the low energy elec- 
tronic density of states. In this Letter, we investigate this 
problem using a suitable generalization of the Momentum 
Average (MA) approximation which allows us to in- 
vestigate non-parabolic bands (tight-binding models) for 
any el-ph coupling strength. This is because this method 
is accurate for any coupling strength, being exact in the 
asymptotic limits of both very weak and very strong cou- 
pling. We demonstrate that the conclusions of Refs. do 
not apply in the low density limit for a tight-binding dis- 
persion. In fact, the exact opposite holds, namely the 
effective el-ph coupling is suppressed by SO coupling. 



We also calculate spin-dependent spectral functions and 
show that in a certain regime, the coherent polaron band 
is dominated by contributions from only one SO elec- 
tronic band ( the "-" band). This has major implications 
for spin dependent transport through such materials, al- 
lowing for the possibility to manipulate the effective mass 
and the spin polarization of quasiparticles by tuning the 
el-ph and SO interactions. 

We consider a single electron on a two dimensional 
square lattice with SO coupling, which also interacts with 
optical phonons of energy U {h ~ 1) through a local 
Holstein-type interaction [lo| . The Hamiltonian can be 
written in terms of k-space spinors as: 
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where = — 2t[cos(fc.^a) + cos(fcj,a)] is the 2D free- 
electron dispersion for nearest-neighbor hopping and 
(/)k = 2Vs[i sin(A:a;a) -I- sin(fcya)] describes the Rashba SO 
coupling. Different dispersions and/or SO couplings can 
be studied similarly. As usual, c]^ ^ is the creation oper- 
ator for an electron with momentum k and spin cr, while 

creates a phonon of momentum q. Both electron spin 
channels interact with the phonons through the local lat- 
tice displacement, with an el-ph coupling constant g. The 
lattice has a total of TV sites and periodic boundary con- 
ditions, and in the end we let N — > 00. Momentum sums 
arc over the Brillouin zone (BZ). 

The 2x2 Green's functions for the non-interacting 
system (no el-ph coupling) are defined as: 



Go(k,c^) 



where Go{uj) = [lu 
responding to Tip = 
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Go+(k,u.) e'^-Go-{k,uj) 
'«-Go-(k,c^) Go+(k,c^) 



(1) 



- Ho + iv] ^ is the resolvent cor- 
7i|g=o- The symmetric and anti- 
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symmetric Green's functions are: 

Go± (k, w) = - ( 7-— 7 ± j—— 7 I 

2 + 17] - Ek + \(Pk\ UJ+ 11] - Cfe - |0fc| / 

and the phase factor is = 0k/|'/'k|- To avoid confusion 
with scalars, all 2 x 2 matrices will be identified by a bar, 
hence the Go(k, a;) notation. 

The full Green's function, defined as usually: 

G{k,u;)^{0\(^'^Giu;)(cl 4,) |0), (2) 

can in principle be calculated exactly by applying the 
equation of motion method. The full Green's function 
can be related to the non-interacting one by using the 
Dyson equation, G{lu) = Go(a;) + G{uj)VGo{uj) - here 
V — H — 'Ho is the el-ph interaction. As shown previously 
for the Holstein Hamiltonian [ll|, [l^] , the repeated use 
of the Dyson identity generates a system of equations in- 
volving generalized Green's functions with various num- 
bers of phonons. In the presence of SO interactions, the 
equations of motion are similar, except they are in terms 
of 2 X 2 Green's functions. As a result, we can implement 
the MA approximations in the same way, and all conclu- 
sions regarding accuracy for all el-ph coupling strengths, 
sum rules obeyed exactly by the spectral weight, etc. re- 
main true. A somewhat analogous procedure was used 
to compute the 2x2 Green's function describing rippled 
graphene however there the matrices are related to 
different sublattices, not to different spin projections. 

For all levels of the MA approximation, the Green's 
function can be written in the standard form: 

G(k, uj) = [Go(k, uj)-^ - S(k, Lu)]-\ (3) 

where the self-energy I](k, w) has different expressions 
depending on the level of MA approximation used. For 
the simplest, least accurate MA'"^ level [9, the self- 
energy has no k dependence and is given by an infinite 
continued fraction Ejy,j^(o) (w) = defined by 

An{uj) = nga{uj - nQ,)[l - g^go{uj - nn)A„+iiuj)]~^, (4) 

where 

.-.H^EGolk,.,. » ) (5) 

k ^ ^ 

is the momentum average of the noninteracting Green's 
function over the BZ. Note that because the off-diagonal 
part is anti-symmetric, its average over the Brillouin zone 
vanishes, thus go{i^) and all An{uj) matrices are diagonal. 

As discussed extensively in Refs. [ll| and MA^°^ is 
accurate for ground state properties but it fails to prop- 
erly predict the polaron-|-one phonon continuum. As a 
result, it overestimates the polaron bandwith. This prob- 
lem is fixed at the MA^-'^' level, where a phonon is allowed 



to appear away from the polaron cloud. For the Holstein 
model (and by extension, in the presence of SO coupling) 
both these approximations predict k-independent self- 
energies. Here, they are diagonal as well, i.e. phonon 
emission/absorption is not allowed to scatter the elec- 
tron between the two spin-polarized bands. These, of 
course, are approximations, although it is worth point- 
ing that all non-crossed diagrams give k-independent and 
diagonal contributions to the self-energy. This is why 
the self-consistent Born approximation prediction for the 
self-energy is also k-independent and diagonal Q. 

In MA, however, the effect of non-crossed diagrams 
is included for MA^^^ and higher levels (in variational 
terms, these allow two or more phonons to appear away 
from the main polaron cloud and the order of their cre- 
ation/annihilation now becomes relevant). We therefore 
report MA*^^^ results here. Followin g a similar procedure 
to that described in detail in Ref. [I4I, the MA^^^ self- 
energy is found as: 

^MA(2)(k,Ci^) 5(0), (6) 

given by the solution of the system of coupled equations 
for the unknown 2x2 matrices x{i): 

Y,M.A^,Lu)xU) = e^'^^^g^Goi-i^^). (7) 
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The sum is over lattice sites i = (ix, iy) located at = 
ixix + iyay. The 2x2 matrices Mij(\i,Lu) are: 

Afoo = 1 - g'go(w)5o('i) {2d^i - a^i) , (8) 
M,o = -~g^goi^)e''' '^'GoH, ^) {2a^,' - a^,') (9) 

for i 7^ 0, and for both i,j 7^ 0: 

A~U, = a2i5.,,l - g^e'^^-^^GoU,^) 

X [(^2 - Ai)(5.,-, + Goi-i - j,i^)a^i] ■ (10) 

Here we defined = 1 — g^go{Lij){Ai — Aj) where Ai = 
Ai{lu — 2n), A2 = A2{llj — il) , A3 = ^3(0;) are continuos 
fractions defined by Eq. (|4]), and a) = lj — 2ri — 1), 
a) = w - g^go(w)|(i,i)(a^/)|(i4) (since the A, a, go matri- 
ces are proportional to 1, the (2,2) diagonal matrix ele- 
ment can be used just as well in the definitions of uj,oj). 
Finally, the real space Green's functions appearing in the 
inhomogeneous terms are given, as usual, by: 

Go(z,c.) = l^e'''^'Go(k,.;). (11) 

k 

It is important to note that for i ^ 0, Go(i,a;) acquires 
off-diagonal components, which lead to the off-diagonal 
contributions in Sjy^^(2) (k, w). Because below the free- 
electron continuum the Go(i,<x') decrease exponentially 
as |Ri| increases, the system in Eq. ([7]) can be truncated 
at a small We truncate at |Ri| ~ 10a, such that the 
relative error of the spectral function is less that 10""^. 
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FIG. 1: (color online) The ground state energy as a function 
of Rashba SO coupling for three values of the el-ph coupling. 
The inset shows the effective mass on a logarithmic scale as 
a function of Rashba SO coupling. 



that effectively decreases. As a result, away from the 
hmit where both A and Vg/t are very small, an increase 
in the SO coupling leads to a drop in the effective mass, 
making it possible to tune the mass of the polaron be- 
tween heavy and light. Light polarons are thus found for 
either small A irrespective of Vg/t, or at large A and large 
enough Vg/t. As we show next, however, their nature 
and spin-character may be very different. 

The MA(2) approximation is quantitatively accurate 
not only for the ground state, but also for high energy 
states, as it satisfies exactly the first 8 spectral weight 
sum rules and with good accuracy all other ones. It is 
thus possible to have an accurate depiction of the high en- 
ergy states by calculating the spectral function. Because 
we expect to have spin polarized bands it is necessary to 
calculate a spin dependent spectral function: 
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/m(Tr[aG(k,a;)]), 



(12) 



Once the self-energy is known, the full Green's func- 
tion can be calculated and in turn provide accurate esti- 
mates for spectral weights, the ground state energy, the 
effective mass etc. In the non-interacting case {g = 0) 
the ground state consists of four degenerate points in k- 
space, (±fcmi„, zbfcmm), where krama = arctan[ys/ {V^t)], 
with the ground-state energy = — 4t cos(fcmi„a) — 
VsV8| sin(fc„ij,ia)|. On the other hand, in the absence 
of SO coupling {Vs = 0), as the effective el-ph coupling 
A = /{AtVL) is turned on, there is a crossover from 
a light, large-polaron to a very heavy, small-polaron at 
A ~ 1. In Fig. [1] we show the ground state energy mea- 
sured from for weak, medium and strong effective el- 
ph couplings, as a function of the Rashba SO coupling. 
For large SO coupling, the renormalization of the en- 
ergy and effective mass (shown in the inset) is strongly 
suppressed, indicating large-polaron behavior even when 
A = 2. This contradicts reported results for a parabolic 
band dispersion 0], which are based on the fact that 
the density of states at the band edge has a square root 
singularity, because in a continuum model the locus of 
momenta defining the ground state is a circle of radius 
kmin ■ This is not gcnerically true for a tight binding dis- 
persion, where the Van Hove singularity is shifted from 
the 4-point degenerate ground state to higher energies. 
As expected, though, the results for a parabolic disper- 
sion can be recovered by our tight-binding model in the 
regime where both A and Vs/t are very small. Indeed, 
for A = 0.75, the effective mass increases slightly with Vs 
at small K,, before decreasing at larger K, values. Such 
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behavior is more apparent as A ^ 

Our results can be understood by noting that the free- 
electron bandwidth increases with increasing Vg- This 
results in an effective el-ph coupling, which compares the 
polaron binding energy to this renormalized bandwidth, 



where S are the Pauli matrices. The direction of A{k, lo) 
gives the direction of the expectation value of the spin, 
while its magnitude gives the density of states with mo- 
mentum k. We know that in the non-interacting case as 
we go around the F-point, eigcnstatcs have spin perpen- 
dicular to their momentum direction and rotating clock- 
wise for one band and anti-clockwise for the other. For 
the coupled system we observe similar spin eigenstates 
and thus choose to plot the following quantity: 



i(k,Lj) ^ [uk X A{\<i,uj)] 



(13) 



where Uk and Uz are unit vectors parallel to k, respec- 
tively z-axis. The two spin polarized bands will now cor- 
respond to opposite signs of A(k, w). 

Since the polaron bandwidth cannot exceed f2, we ex- 
pect the character of the polaron band to depend on the 
relation between $7 and the energy difference between the 
spin-split electron bands. In order to exemplify this we 
plot i(k,u;) in Fig. [H for A = 1, f7 = QM and = 0.4<. 
The spectral function is shown along the (0, 0)-(7r, tt) line 
in order to intersect the ground state at {kmin, kmin)- In 
this case J7/(£'o~4t) > 1 and we see two coherent polaron 
bands corresponding to the two spin polarizations, with 
similar quasiparticle weights. We conclude that when 
is larger than the SO splitting, the polaronic quasipar- 
ticles are rather similar to the non-interacting electrons, 
except for the renormalized mass and supressed quasi- 
particle weight. Of course, the spectrum above Eqs + ^ 
becomes incoherent due to el-ph scaterring. 

In Fig. [3] we plot i(k,cj) for VL = 0.2t, A = 1.0 and 
Vs = QM. Now Vl/{EQ~At) < 1 and because the polaron 
bandwidth cannot exceed O, it is dominated by the "-" 
band. There is a large difference between the quasiparti- 
cle weights of the two coherent polaron bands (note the 
different positive and negative scales for the the contour 
plot), and the effective mass of the dominant "-" band 
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FIG. 2: (color online) Spin dependent spectral function 
A(k,w), for A = 1, n = 0.8t and Vs = OAt. The right panel 
shows only the "+" band for clarity. 



is much smaller than that of the low-weight "+" band. 
Higher energy states have small weights and are highly 
incoherent, i.e. short lived. 

Consider now injection of a current in such a system. 
Whereas in a regime like that depicted in Fig. [2]we expect 
the spin to precess between the two coherent polaronic 
bands, like it does for non-interacting electrons 1^, lit , 
in a regime like in Fig. [3] only the spin-component parallel 
to the "-" band can be efficiently transmitted through 
the system, which therefore acts as an intrinsic "spin- 
polarizer" . 

This becomes more and more apparent as one moves 
further into the asymptotic limit where both the SO 
and the el-ph couplings are large compared to il, i.e. 
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FIG. 3: (color online) Spin dependent spectral function 
l(k,a;), for A = 1, SI = 0.2t and Vs = 0.8t. The right-hand 
side panel shows only the "-f" band for clarity. 



{Eo - 4t) > and A > 1. In this limit the "-" band in 
the coherent polaron spectrum becomes lighter and has 
a higher quasiparticle weight, whereas the band es- 
sentially vanishes from the coherent spectrum (its quasi- 
particle weight is extremely low and its effective mass is 
extremely large) . The resulting light polaron is thus very 
different from the one we find in the small A regime, which 
has both bands in the coherent spectrum with roughly 
equal quasiparticle weights and effective mass. 

This demonstrates that an interplay between SO and 
el-ph couplings allows indeed for different tailoring of the 
properties of the two spin-polarized bands, whereas one 
is well described by a long-lived, light quasiparticle while 
the other becomes highly incoherent. This will natu- 
rally lead to very different conductivities for the two spin 
polarizations, making such a material ideal as a source 
and/or detector of spin-polarized currents - these are im- 
portant components needed for many spintronics appli- 
cations. These conclusions arc based on the use of the 
MA approximation which is known to be highly accurate 
for all el-ph couplings. Moreover, at the MA'^' level we 
use here, it results in a non-diagonal, k-dependent self- 
energy, therefore our results properly include phonon- 
mediated scattering between the two electronic bands. 
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